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Abstract 

We present the number of dimers Nd{n) on the Sierpinski gasket SGd{n) at stage n with di- 
mension d equal to two, three, four or five, where one of the outmost vertices is not covered when 
the number of vertices v{n) is an odd number. The entropy of absorption of diatomic molecules 
per site, defined as SsGd — li™n-»oo In A''rf(n)/t;(n), is calculated to be ln(2)/3 exactly for SG2{n). 
The numbers of dimers on the generalized Sierpinski gasket SGd,b{n) with d = 2 and 6 = 3, 4, 5 
are also obtained exactly. Their entropies are equal to ln(6)/7, ln(28)/12, ln(200)/18, respectively. 
The upper and lower bounds for the entropy are derived in terms of the results at a certain stage 
for SGd{n) with d = 3,4,5. As the difference between these bounds converges quickly to zero as 
the calculated stage increases, the numerical value of SsG^i with d = 3, 4, 5 can be evaluated with 
more than a hundred significant figures accurate. 
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I. INTRODUCTION 



The enumeration of close-packed dimers N{G) on a graph G was first considered by 

bwler and Rushbrooke in enumerating the absorption of diatomic molecules on a surface 
1 

11]. The dimer coverings of a graph is a classical model in statistical physics and is called 
perfect matchings in mathematical literature. The dimer model on the square lattice was 
solved exactly by Kasteleyn [2] and Temperley and Fisher 0, 4|. The model is equivalent to 
various other statistical mechanical problems. For example, the zero-field partition function 
of Ising model on a planar lattice can be formulated as a dimer model on an associated 
planar lattice It is also well known that there is a bijection between close-packed 

dimer coverings and spanning tree configurations on two related planar lattices j3] • A recent 
review on the enumeration of close-packed dimers on two-dimensional regular lattices is 
summarized in Ref. 8j. It is of interest to consider dimer coverings on self-similar fractal 
lattices which have scaling invariance rather than translational invariance. Fractals are 
geometric structures of noninteger Hausdorff dimension realized by repeated construction of 
an elementary shape on progressively smaller length scales jol, [lo|. A well-known example 



of frac 
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tal is the Sierpinski gasket w 
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lich has been extensively studied in several contexts 
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2l[ |. A dimer coverings will leave at least one vertex 



uncovered when the total number o: 
with both length and width odd 22 



vertices is an odd number, e.g., the rectangular lattice 



231]. The vacancies that are not covered by any dimers 



can be considered as occupied by monomers. We allow such possible vacancies occur on the 
outmost vertices of the Sierpinski gasket. We shall derive rigorously the numbers of dimer 
coverings on the two-dimensional Sierpinski gasket and its generalization, and obtain upper 
and lower bounds for the entropy on the Sierpinski gasket with dimension equal to three, 
four or five. 



II. PRELIMINARIES 

We first recall some relevant definitions in this section. A connected graph (without 
loops) G = (y, E) is defined by its vertex (site) and edge (bond) sets V and E 



24 



25]. 



Let v{G) = \V\ be the number of vertices and e{G) = \E\ the number of edges in G. The 
degree or coordination number ki of a vertex Vi E V is the number of edges attached to it. 
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A fc-regular graph is a graph with the property that each of its vertices has the same degree 
k. In general, one can associate a dimer (monomer) weight to each dimer (monomer) (see, 



for example 22|). For simplicity, all dimer (monomer) weights are set to one throughout 
this paper. 

When the size of the graph increases as v{G) oo, the number of dimer coverings N{G) 
grows exponentially in v{G). The entropy of absorption of diatomic molecules per site is 
given by 

IniV(G) 

Sg = lim — , 2.1 

v{G)~*oo v{G) 

where G, when used as a subscript in this manner, implicitly refers to the thermodynamic 
limit. The dimer coverings considered here may not be close-packed dimers since there may 
be vacancies on the outmost vertices as mentioned above. Notice that we define the entropy 
per site rather than entropy per dimer. They differ by a factor of two in the thermodynamic 
limit regardless the presence of vacancies on the outmost vertices. 

The construction of the two-dimensional Sierpinski gasket SG2{n) at stage n is shown 
in Fig. [H At stage n = 0, it is an equilateral triangle; while stage n + 1 is obtained by 
the juxtaposition of three ra-stage structures. In general, the Sierpinski gaskets SGd can be 



built in any Euclidean dimension d with fractal dimensionality D = \n{d+ l)/ln2 12|. For 
the Sierpinski gasket SGd{n), the numbers of edges and vertices are given by 

e{SG,{n)) = ("^ 2 ^ = 1^"^ + ' ^^'^^ 

v{SG,{n)) = ^^[{d + ir + l] . (2.3) 

Except the {d + 1) outmost vertices which have degree d, all other vertices of SGd{n) have 
degree 2d. In the large n limit, SGd is 2(i-regular. 



A 





SG2(0) SG2{1) SG2{2) 5G2(3) 
FIG. 1: The first four stages n = 0, 1, 2, 3 of the two-dimensional Sierpinski gasket SG2{n). 
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The Sierpinski gasket can be generalized, denoted as SGd^b{n), by introducing the side 
length b which is an integer larger or equal to two [26] . The generalized Sierpinski gasket at 
stage + 1 is constructed with b layers of stage n hypertetrahedrons. The two-dimensional 
SG2,bi'n) with 6 = 3 at stage n = 1,2 are illustrated in Fig. [21 and those with 6 = 4, 5 at 
stage n = 1 in Fig. [3l The ordinary Sierpinski gasket SGd{n) corresponds to the 6 = 2 case, 
where the index b is neglected for simplicity. The Hausdorff dimension for SGd.b is given by 
D = \n ln6 26]. For the two-dimensional Sierpinski gasket SG2,biji) that will be 

considered here, the numbers of edges and vertices are given by 

-6(6 + 1) 



e{SG2M) 



v{SG2An)) 



6 + 4 
6 + 2 



6(6+1)]" , 2(6 + 1) 
6 + 2 



+ 



Notice that SGd.b is not fc-regular even in the thermodynamic limit. 



(2.4) 
(2.5) 




SG2dl 



SG2,z{2) 



FIG. 2: The generalized two-dimensional Sierpinski gasket SG2,b{n) with 6 = 3 at stage n — 1^2. 




5G2,4(1) SG2,5(1 



FIG. 3: The generalized two-dimensional Sierpinski gasket SG2,b{n) with 6 = 4, 5 at stage n — 1. 
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III. THE NUMBER OF DIMER COVERINGS ON SG2,b{n) WITH 6 = 2,3,4,5 

In this section we derive rigorously the numbers of dimer coverings on the two-dimensional 
Sierpinski gasket SG2{n), equivalently SG2,2{n), and the generalized SG2,b{n) with b = 
3, 4, 5. Let us start with the definitions of the quantities to be used. They are illustrated in 
Fig. m where only the outmost vertices of SG2fiin) are shown. 

Definition III.l Consider the generalized two-dimensional Sierpinski gasket SG2,b{n) at 
stage n. (i) Define f2,b{n) as the number of dimer coverings such that the three outmost 
vertices are vacant, (ii) Define g2,b{i^) the numbers of dimer coverings such that one 
certain outmost vertex, say the topmost vertex as illustrated in Fig. ^ is occupied by a 
dimer while the other two outmost vertices are vacant. (Hi) Define h2,b{n) as the numbers of 
dimer coverings such that one certain outmost vertex, say the topmost vertex as illustrated 
in Fig. is vacant while the other two outmost vertices are occupied by dimers. (iv) Define 
t2,b{n) as the number of dimer coverings such that all three outmost vertices are occupied by 
dimers. 



FIG. 4: Illustration for the configurations f2.bin), g2.b{n), h2.bin), and t2.b{n). Only the three outmost 
vertices are shown explicitly, where each open circle is vacant and each solid circle is occupied by a dimer. 

A. SG2{n) 

For the ordinary two-dimensional Sierpinski gasket, we use the notations f2{n), g2{n), 
h2{n), and t2{n) for simplicity. Because of rotational symmetry, there are three possible 
g2{n) and three possible h2{n) for non- negative integer n. The initial values at stage zero 
are /2(0) = 1, 5'2(0) = 0, /i2(0) = 1, ^2(0) = 0. The values at stage one are /2(1) = 0, 
(?2(1) = 2, /i2(l) = 0, ^2(1) = 2. The value zero indicates that no such configurations are 
allowed. By Eq. fl2.3p . we have 




f2,b{n) 



92,b{n) 



h2,b{n) 



t2,b{n) 
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j=i VJ/ vj/ 

where the Binomial expansion is used for 3" = (2 + 1)", such that the number of vertices 
for SG2{n) is odd for even n and even for odd n. Therefore, f{n), h{n) are always zero for 
odd n and g{n), t{n) are always zero for even n. Let us denote odd n as 2m + 1 and even n 
as 2m with non-negative integer m in the following discussion for SG2{n). These quantities 
satisfy simple recursion relations. 

Lemma III.l For any odd n = 2m + 1 > 0, 

/2(2m + 2) = 2^^^2m + l) , (3.2) 

/i2(2m + 2) = 2^^(2m + l)t2(2m + 1) . (3.3) 

For any even n = 2m > 0, 

g2{2m + 1) = 2/2 (2m) /i^ (2m) , (3.4) 

t2(2m + 1) = 2hl{2m) . (3.5) 

Proof The Sierpinski gasket SG2{n + 1) is composed of three SG2{n) with three pairs 
of vertices identified. For each pair of identified vertices, either one of them is originally 
occupied by a dimer while the other one is vacant. The number /2(2m + 2) for SG2{2m + 2) 
consists of two configurations where all three of the S'G2(2m + 1) are in the 5'2(2m + 1) status 
as illustrated in Fig. such that Eq. fl3.2p is verified. 




FIG. 5: Illustration for the expression of /2(2m + 2). 

Similarly, /i2(2m + 2) and 5'2(2m + 1), t2{2m + 1) can be obtained with appropriate 
configurations of its three constituting blocks as illustrated in Figs. El [7] and [H] to verify Eqs. 
f lOl) . M and ([33]), respectively. □ 
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FIG. 6: Illustration for the expression of ft,2(2m + 2). 




FIG. 7: Illustration for the expression of 52 (2m + 1). 
It is elementary to solve f2{n), g2{n), h2{n), t2{n) in order to obtain the entropy for 5^2- 



Theorem III. 2 For the two-dimensional Sierpinski gasket SG2{n) at stage n = 2m or 
n = 2m + 1, 

f2(2m) = h2(2m) = 2^2(2^ 

(3.6) 

/2(2m + 1) = /i2(2m + 1) = 

g2{2m) = t2(2m) = 
^2(2m + 1) = t2(2m + 1) = 2^2{2m+i) 
where the exponent is 

72(n) = ^(3"-l) . (3.8) 

Define the number of dimer coverings N{SG2{n)) in Eq. h2. 1\) equal to f2{n = 2m) and 
equal to g2{n = 2m + 1) for even and odd n, respectively. With v{SG2{n)) = |(3" + 1), the 
entropy is given by 

SsG2 = -ln2 ~ 0.23104906018... (3.9) 
3 




FIG. 8: Illustration for the expression of t2{'2m + 1). 
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B. SG2,3{n) 

For the generalized two-dimensional Sierpinski gasket 5'G2,b(^) with 6 = 3, we have 
viSG^M) = li^Lll = e'^ + 2 + ^(6" - 1) = 6" + 2 + 2 Q 5^~' (3.10) 

by Eq. (12.51) . such that the number of vertices is equal to three for n = and becomes even 
for all positive integer n. Therefore, f2,3{n) and /i2,3(^) are always zero for positive integer 
n, while the initial values remain /2,3(0) = 1, (72,3(0) = 0, /i2,3(0) = 1 and ^2,3(0) = 0. g2,z{ji) 
and t2^3{n) satisfy recursion relations. 

Lemma III. 3 For any positive integer n, 

g2,3{n + l) = Qgl^{n)t2,3{n) , (3.11) 
t2,3{n + l) = Qgl^{n)tl^{n) , (3.12) 

and for n = 0, 

^?2,3(1) = 6/|3(0)/i^,3(0) = 6 , (3.13) 

^2,3(1) = 6/2,3(0)/l^,3(0) = 6 . (3.14) 

Proof The Sierpinski gasket SG2,3{n+ 1) is composed of six SG2,z{n) with six pairs of 
vertices identified and a set of three vertices identified. For each pair of identified vertices, 
either one of them is originally occupied by a dimer while the other one is vacant. For the 
set of three vertices, one of them is originally occupied by a dimer while the other two are 
vacant. The number g2,3{n + 1) for positive n consists of six configurations where five of the 
SG2,'i{n) are in the (72,3 (^) status and one in the t2,3(^) status as illustrated in Fig. [9|, such 
that Eq. (I3.1ip is verified. 




FIG. 9: Illustration for the expression of 52,3 (" + 1) with positive n. The multiplication of two on the 
right-hand-side corresponds to the reflection symmetry with respect to the central vertical axis. 
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Similarly, t2,3{n + 1) with positive n for SG2,z{n + 1) can be obtained with appropri- 
ate configurations of its six constituting 5* 6*2,3 1'^) as illustrated in Fig. [10] to verify Eq. 
fl3:T2D . Finally, c/2,3(l) and ^2,3(1) in Eqs. fl313l) and fl3:T4ll are verified by Figs. [H] and [121 
respectively. □ 



FIG. 10: Illustration for the expression of t23{n + 1) with positive n. The multiplication of three on the 
right-hand-side corresponds to the three possible orientations of SG2,3{n, + 1) 




FIG. 11: Illustration for the expression of 52,3(1)- The multiplication of two on the right-hand-side corre- 
sponds to the reflection symmetry with respect to the central vertical axis. 



FIG. 12: Illustration for the expression of ^2,3(1)- The multiplication of three on the right-hand-side 
corresponds to the three possible orientations of SG2,3{n + 1) 

It is elementary to solve g2,3{n) and t2^3{n) for positive n in order to obtain the entropy 
for SG2,3- 

Theorem III. 4 For the generalized two-dimensional Sierpinski gasket SG2,3{n) at stage 
n> 0, 





^2,3(n)=t2,3(n) = 6^^-^('^) , 



(3.15) 
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where the exponent is 

72,3H = J(6"-l) . (3.16) 

Define the number of dimer coverings N{SG2,3{n)) in Eq. ^2. 1\) equal to t2.3{n). With 
v{SG2,3{n)) = (7(6)" + 8)/5, the entropy is given by 

Ssa^ .^ = ^ln6 ~ 0.25596563846... (3.17) 



C. 5G2,4(n) 

For the generalized two-dimensional Sierpinski gasket SG2fi{n) with 6 = 4, we have 

,.^SGUn)) = '-^^ = 3 + 1(10" - 1) = 3 + 1 1 ('> (3.18) 

j=i V-^/ 

by Eq. (12.51) . such that the number of vertices is always odd for any n. Therefore, g2,i{n) and 
^2,4 (^) are zero for all n, while the initial values remain /2,4(0) = 1, ^'2,4(0) = 0, /i2,4(0) = 1 
and ^2,4(0) = 0. f2,i{n) and /i2,4(^) satisfy recursion relations. 

Lemma III. 5 For any non-negative integer n, 

f2A{n + l) = 28fl,{n)hl,{n) , (3.19) 

/i2 4(n + 1) = 28fl^in)hl^{n) . (3.20) 

Proof The Sierpinski gasket SG2,A{n + 1) is composed of ten S'G'2,4('^) with nine pairs 
of vertices identified and three sets of three vertices identified. The number /2,4(^ + 1) for 
non-negative n consists of twenty eight configurations where four of the 5'G'2,4(^) are in the 
/2,4(n) status and six in the h2,i{n) status as illustrated in Fig. [131 such that Eq. (13.191) is 
verified. 

Similarly, /i2,4(^+l) with non-negative n for SG2,i{n+l) can be obtained with appropriate 
configurations of its ten constituting S'G2,4('^) as illustrated in Fig. [Tl]to verify Eq. (13.201) . 

□ 

It is elementary to solve f2,i{n) and h2,i{n) in order to obtain the entropy for S'G'2,4- 
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FIG. 13: Illustration for the expression of f2A{n + 1) with non-negative n. The multiplication of three on 
the right-hand-side corresponds to the three possible orientations of SG2,4{n + 1). 

Theorem III. 6 For the generalized two-dimensional Sierpinski gasket SG2,4,{n) with non- 
negative integer n, 

f2A{n) = h2,A{n) = 28^^-*(") , (3.21) 

where the exponent is 

72,4H = 1(10^-1) . (3.22) 

Define the number of dimer coverings N{SG2,i{n)) in Eq. Ii2. 1\) equal to f2,i{n). With 
v{SG2,'i{n)) = (4(10)" + 5)/3, the entropy is given by 

SsG2A = — ln28 ~ 0.27768370918... (3.23) 



D. 5G2,5(n) 

For the generahzed two-dimensional Sierpinski gasket SG2,b{n) with 6 = 5, we have 
v{SG2,,{n)) = ^^^^^^^^^ = 15« + 2 + ^(15" - 1) = 15" + 2 + ^ f'') W (3.24) 
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FIG. 14: Illustration for the expression of h2,4:{n + 1) with non-negative n. The multiplication of two on 
the right-hand-side corresponds to the reflection symmetry with respect to the central vertical axis. 

by Eq. fl2.5l) . such that the number of vertices is always odd for any n. Therefore, 5'2,5(^) and 
t2,5{n) are zero for all n, while the initial values remain /2,5(0) = 1, fi'2,5(0) = 0, /i2,5(0) = 1 
and ^2,5(0) = 0. The figures of the recursion relations for /2,5(^) and h2^^{n) are too many 
to be shown here, and we state the following Lemma without proof. 

Lemma III. 7 For any non-negative integer n, 

/2,5(n + 1) = '2mfl^{n)hl^{n) , (3.25) 
/i2,5(n + 1) = 200/2%(n)/i^°5(n) . (3.26) 

It is elementary to solve f2,5{n) and /i2,5('^) in order to obtain the entropy for 5^2,5 • 
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Theorem III. 8 For the generalized two-dimensional Sierpinski gasket SG2^5{n) with non- 
negative integer n, 

/2,5H = /i2,5H= 200^^'^(") , (3.27) 

where the exponent is 

72,5(n) = ^(15" - 1) . (3.28) 

Define the number of dimer coverings N{SG2,5{n)) in Eq. \2. 1\) equal to f2,5{n). With 
v{SG2,5{n)) = (9(15)" + 12)/7, the entropy is given by 

SsG^, = ^ In 200 ~ 0.29435096480... (3.29) 
18 

As the generahzed two-dimensional Sierpinski gasket SG2,b{n) for any h is planar, it 
appears that the number of dimer coverings can be solved exactly. However, the number of 
configurations to be considered increases as h increases and the recursion relations must be 
derived individually for each h. We have been unable to obtain a general expression of the 
number of dimer coverings on SG2,b{n) for arbitrary b. 

IV. THE NUMBER OF DIMER COVERINGS ON SGd{n) WITH d = 3,4,5 

In this section we present the number of dimer coverings on the Sierpinski gasket SGd{n) 
with d = 3,4,5 which is not planar. Instead of solving exactly the entropies for these 
Sierpinski gaskets, we obtain accurate upper and lower bounds for them. 

A. 5G3(n) 

For the three-dimensional Sierpinski gasket SGs{n), we use the following definitions. 

Definition IV. 1 Consider the three-dimensional Sierpinski gasket SG^ln) at stage n. (i) 
Define fzi^) as the number of dimer coverings such that the four outmost vertices are vacant, 
(a) Define h-i{n) as the number of dimer coverings such that two certain outmost vertices 
are occupied by dimers and the other two outmost vertices are vacant. (Hi) Define S3(n) as 
the number of dimer coverings such that all four outmost vertices are occupied by dimers. 
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As the number of vertices for SG^ln) is always even by Eq. (12. 3p . we do not have the 
dimer coverings such that one certain outmost vertices is occupied by a dimer and the other 
three outmost vertices are vacant, or one certain outmost vertices is vacant and the other 
three outmost vertices are occupied by dimers. The quantities f^iji), h^{n), and s^in) are 
illustrated in Fig. [151 where only the outmost vertices are shown. There are (2) = 6 
equivalent h^{n). The initial values at stage zero are /3(0) = 1, h^{Q) = 1, ^3(0) = 3. These 
quantities satisfy recursion relations. 




fsin) h^{n) S2,{n) 

FIG. 15: Illustration for the dimer coverings fz{n)^ h^{n) and 53(71). Only the four outmost vertices are 
shown explicitly, where each open circle is vacant and each solid circle is occupied by a dimer. 



Lemma IV. 1 For any non-negative integer n, 

h{n + l) = Sh{n)hl{n) , (4.1) 

h^{n + 1) = AUn)hl{n)s;{n) + Ahl{n) , (4.2) 
S3(n + 1) = 8/i^(n)s3(n) . (4.3) 

Proof The Sierpinski gasket SG^{n + 1) is composed of four SG^{n) with six pairs 
of vertices identified. For each pair of identified vertices, either one of them is originally 
occupied by a dimer while the other one is vacant. The number /3(r;, + 1) for non-negative 
n consists of eight configurations where one of the SG^{n) are in the f^in) status and the 
other three are in the h^{n) status as illustrated in Fig. [161 such that Eq. (14. ip is verified. 



Similarly, h^{n + 1) and 53(71, + 1) for SG^{n + 1) can be obtained with appropriate 
configurations of its four constituting SG^{n) as illustrated in Figs. [T71 and [18] to verify 
Eqs. (53) and (jlJD, respectively. □ 



Dimer coverings on the Sierpinski gasket 



16 





x4 + 




FIG. 16: Illustration for the expression of f3{n + 1). The multiplication of four on the right-hand-side 
corresponds to the four possible orientations of SGsin + 1). 




+ 




x2 + 



x2 + 




FIG. 17: Illustration for the expression of /i3(n + 1). The multiplication of two on the right-hand-side 
corresponds to the reflection symmetry with respect to the central vertical axis. 



The values of fsin), hsln), 53(72) for small n can be evaluated recursively by Eqs. (14. ip - 
f l4.3p . but they grow exponentially, and do not have simple integer factorizations. To estimate 
the value of entropy for SG3, we define the ratio 



and its limit 



^3 = lim 03 (n) . 



(4.4) 



(4.5) 



Lemma IV. 2 Sequence {a3(?T,)}^^ decreases monotonically. The limit is equal to v^. 




FIG. 18: Illustration for the expression of s^{n + 1). The multiplication of four on the right-hand-side 
corresponds to the four possible orientations of S 03(11 + 1). 
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Proof From Eqs. (14. ip and (14.31) . the ratio 33(71) / f^ln) is invariant, that is equal to 
S3 (0)7/3(0) = 3. Eq. (H^D can be modified to be 

h,in + 1) = 12fi{n)hlin) + 4/i^(n) . (4.6) 

Although 03(0) = 1, it is clear that as{n) is bounded below by for positive integer n 
because 

hlin + 1) - 3/|(n + 1) = [12flin)hlin) - Ahtin)]' > . (4.7) 
It follows that 0:3 (n) decreases for positive n because 

hsjn) _ hjn + l) ^ hlin) - 3/|(n) ^ 

Mn) Mn + 1) 2f,{n)h{n) " ' ^ ' ^ 

which implies that the limit as exists. From Eqs. (14. ip and (14. 6p . we have 

hjn + l) _3Mn) Ihjn) 
h{n + l) 2h{n) 2h{n)- 

By taking the large n limit in Eq. (14. 9p . 0:3 is solved to be -\/3. □ 

The general expressions for fz{n) and ^{n) in terms of quantities at stage m < n can be 
written as follows. 

Lemma IV. 3 For a non-negative integer m and any positive integer n > m, 

2(4)"-"'+l-5-3(-l)"~'" 2(4)"-™+3(-l)"~"' 3(4-|n-m _3(_^)n-m 

fsin) = 2 10 /3(m) 5 haim) 5 

n—m ^ ^ 3(4)J~^ -3(-i)J~-'- 

5 



n [s + alin-j) 



X 

i=2 



(4.10) 



hsin) = 2 5 f^(m) 5 /i3(m) s 

n—m 3(4)J~^+2(-l)J~^ 

X n [3 + "3(^-j)] ' • (4.11) 

Here when n — m = 1, the product with lower limit two is defined to be one. 
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Proof It is clear that Eqs. fl4.10p and (14. lip are valid for n = m + 1 since fsim + 1) = 
8f3{m)hl{m) and h^^m + l) = 4:f^{m)hl{m)[3 + al{m)] by Eqs. (14.11) and (14. 6p . respectively. 
Consider Eq. (I4.10p holds for a certain positive integer n > m, then 



f3in + l)=8Mn)hl{n) 



2(4)"-'" + l-5-3(-l)"~'" 2(4)"-'"+3(-l)"~'" 3(4)"-m _3j^_j^^n-m 

8x2 10 fsim) ^ hs{m) s 



n—m 
J=2 



— 5^ 3(4)"-'"+l-15+3{-l) 

2 5 



n — m 6(4)^^^"' — 6( — 1)^^"""^ 



x/i3(m) 



9(4)"-'-"-|-6(-l)"-'" 



n—m 

n [3 + a2(n-j) 



_^ 9(4)J"l+6(-l)J~l 
5 



8(4)"-'"+! -5+3(-l)"-"' 8(4)"-'"-3(-l)"-'" 12(4)"-'"+3(-l)"-™ 

2 10 /3(m) 5 /i3(m) 5 



X H \3 + al{n-j) 
i=2 



^ 12(4)J-l+3(-l)J-l 

5 



3 + al{n - 1) 



'him)' 



n—m+l 



X 



n [3 + «^(n + l-j) 



_^ 3(4)J-l-3(-lH-l 
5 



J=2 



(4.12) 



such that Eq. (I4.1UI) is proved by induction. Eq. (14. lip can be established by the same 
procedure. □ 



Let us state the following lemma without proof. 



Lemma IV.4 //X(n + 1) = for non-negative integer n and Xq is known, then 



X{n) 



c 



2"~1 



(4.13) 



From above lemmas, we have the following bounds for the entropy. 
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Lemma IV. 5 The entropy for the number of dimer coverings on SG^ln) is bounded: 

V3e3(m)3 f21n/3(m) + 31n/i3(m) + 51n2 + ln3 V^€3{m)^ 

720(4)™ - 10(4)'^ ^ 40(4)™ / 



< S, 



m 



i2 



< , (4.14) 

- 40(4)™[2v^-e3(m)] 

where m is a positive integer and 63(77.) is defined as 0:3(77.) — \/3. 

Proof As fsin) only differs from Ss{n) by a factor of three, which is insignificant in the 
definition of the entropy, we will substitute N{G) = fs^n) in Eq. (12.11) for 5^3 so that 

Ssc. = lim . (4.15) 



n—*oo 



By Lemma [lV.3t we have 



2(4)"-™ + 3(-l)"-™ , 3(4)n->n_3(_l)n->n 

ln/3(n) = ln/3(m)H ln/i3(m) 

+ ^ ' In2 + A3(n,m) , (4.16) 



where 



3r4V'"i - 3('-l 

As{n, m) = Y^ ^ ' ln[3 + a^(n - j)] , (4.17) 

i=2 

which is bounded as follows. 

By Lemma [IV. 2t we know 63(77) decreases monotonically to zero for positive integer n. 
£3(1) = ^3(1) - ^ = 2- ^/3. It is easy to find, by Eq. (gS]), that 

e,{n + l) = 0^ . (4.18) 

^ ^ 2(v^ + e3(n)) ^ ' 

Therefore, for any integer n> m with ?t7 fixed, we have 

63(72 + m) = ^^'^;„_^ (1 + 0(77)) (4.19) 

by Lemma [IV. 4 1 and Eq. (14.181) . where o{n) is negative here and o{n) — ^ as ?7 ^ 00. Now 
A3(?7, m) in Eq. (14.171) can be rewritten, 

A,{n, m) = J2 ^ ^ HQ + 2v^e3(n - j) + 63(72 - jf] . (4.20) 

i=2 
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Since 63(71) is small for positive n, the logarithmic term can be written as 



ln[6 + 2V3es{n - j) + e^^n - jf] = In 6 + ^e-sin - j] 



18 



(4.21) 



where ^nj G (0, 1) so that 
A3 (n,m) 



g 3(4)^-^-3(-l)^-^ f^^g ^ V3e,{n-j) ^^ UMn ' j) 



i=2 

n—m 



^^[(4).-i_(_l).-i] 



i=2 

n—m 



i=2 



^nje3(^ - j 

18 



i2-, 



(4.22) 



Because the j = n — m term gives the largest contribution for A3 (n,m), it is easy to see 
that 

4— ^v/363(m)(l-^)(l + o(n)) ^ ^ ^ 4»-'»ln6 ,^ , , 
^ <A3(?2, m) {l + o{n)). (4.23) 

5 5 

On the other hand, since "^^Z^ e-^ln — j)4''^-'^[l — S,n,j^z{f^ ~ j)^/18] is less than 

A=0 



Z]r=o^ ^4" ™ -'^e3(m + i), we have 



4"-"^ In 6,^ , 3(4)"-"^,^ , Ve3(m)x2' 

A3(n,m) (l + o(n)) < ^ ^ ' "^-^^ ^ ' ^ ' 



10 



;i-<«)) E (If) 



< 



3(4)'^-'^e3(m)(l + o(n)) 
10[2v^-e3(m)] 



(4.24) 



where we use Eq. (14.191) and the inequality 

n—m— 2 



n—m— 2 



2^ 

X = X 



n—m— 2 , 



X 



(4.25) 



for any < a; < 1. The proof is completed by taking the infinite n limit in Eq. (I4.15p . □ 

The difference between the upper and lower bounds for Ssc-i quickly converges to zero 
as m increases, and we have the following proposition. 



Proposition IV. 1 The entropy for the number of dimer coverings on the three-dimensional 
Sierpinski gasket SG^{n) in the large n limit is SsCi = 0.42896389912.... 
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By Eq. f l4.19p . we know 




,6 



(4.26) 



such that SsGs can be calculated with more than a hundred significant figures accurate when 
m is equal to seven in Eq. (14.141) . It is too lengthy to be included here and is available from 
the authors on request. 

B. SG^in) 

For the four-dimensional Sierpinski gasket SG4^{n), we use the following definitions. 

Definition IV. 2 Consider the four- dimensional Sierpinski gasket SGi{n) at stage n. (i) 
Define f^in) as the number of dimer coverings such that the five outmost vertices are vacant, 
(a) Define hi{n) as the number of dimer coverings such that two certain outmost vertices 
are occupied by dimers and the other three outmost vertices are vacant. (Hi) Define s^{n) as 
the number of dimer coverings such that one certain outmost vertex is vacant and the other 
four outmost vertices are occupied by dimers. 

By Eq. (12. 3p . we have 



such that the number of vertices for SGi{n) is always odd. Therefore, we do not have the 
dimer coverings such that one certain outmost vertices is occupied by a dimer and the other 
four outmost vertices are vacant, or three certain outmost vertices are occupied by dimers 
and the other two outmost vertices are vacant, or all five outmost vertices are occupied by 
dimers. The quantities fiiji), hi{n), and Si{n) are illustrated in Fig. [191 where only the 
outmost vertices are shown. There are (2) = 10 equivalent h/^{n) and (^) = 5 equivalent 
s^in). The initial values at stage zero are again /4(0) = 1, /i4(0) = 1, 54(0) = 3. 



V 



{SG,{n)) = |(5" + 1) = 2(5)" + 3 + ^(5" - 1) = 2(5)" + 3 + ^ f '')4^ (4.27) 




We write a computer program to obtain following recursion relations. 
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S4(n) 

FIG. 19: Illustration for the dimer coverings fA(n), hi{n), 54(71). Only the five outmost vertices are shown 
explicitly, where each open circle is vacant and each solid circle is occupied by a dimer. 

Lemma IV. 6 For any non-negative integer n, 

U{n + 1) = AQU{n)hl{n)si{n) + 2Ahl{n) , (4.28) 

hi{n + 1) = 2Afi{n)hl{n)sl{n) + AOhl{n)si{n) , (4.29) 
S4(n + 1) = 8fi{n)hi{n)sl{n) + 56hl{n)sl{n) . (4.30) 

The values of fA{n), h^ln), S4{n) for small n can be evaluated recursively by Eqs. fl4.28p - 
fl4.30p . but they grow exponentially, and do not have simple integer factorizations. To 
estimate the value of entropy for SG4, we define the ratios 

and their limits 

04 = lim a^^n) , /34 = lim ^^{n) . (4.32) 

n— ►oo 71— >oo 

Lemma IV. 7 Sequence {ai{n)}'^^^ decreases monotonically while sequence {/34(n)}^^ in- 
creases monotonically. The ratio (3/^{n) / ai{n) for positive n increases monotonically to one. 

Proof From Eqs. (l4:28|) to (ICTD . we find 

hl{n + 1) - fA{n + l)si{n + 1) = 2hQh\{n)si{n)[hl{n) - fi{n)si{n)f > , (4.33) 



such that 04(71) > (34{n) for n > 0. It follows that for positive n, 

h^in) h^in + 1) 24hl{n)[hi{n) - fi{n)sl{n)] 



h{n) h{n + l) Mn)U{n + l) 

which shows that 04 (n) decreases, and 

S4(n) _ S4(n + 1) _ lQhl{n)sl{n)[f4{n)s4{n) - hl{n)] 
hiin) h^in + 1) h4{n)h4{n + 1) 



> , (4.34) 



< , (4.35) 
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which shows that (34{n) increases. From Eqs. fl4.28p to (14.301) . we have 



and 



/4(n + 1) 5/94(n) + 3a4{n) 

S4(n + 1) /?4(n)[/?4(n) + 7a4,{n)] 



(4.36) 



(4.37) 



h^^n + 1) 3Pi{n) + 5a4{n) 

which leads to = (3^ by taking the hmit n ^ oo. The actual value of 04 and jS^ cannot 
be obtained by solving these equations. The numerical results give 



a^ = (3^ = 0.850772150002... 



(4.38) 



where more than a hundred significant figures can be evaluated when stage n in Eq. (I4.3ip 
is equal to seven. □ 

The general expressions for hi[n) and Si{n) in terms of quantities at stage m < n can be 
written as follows. 



Lemma IV. 8 For a non-negative integer m and any positive integer n > m, 



a(5)"""'-3 3(5)^ 

n4[n) = 2 4 ti4[m) 



n—m 

m 

2=1 



5 + 3 



Pijn - i) 
aAn — i) 



3(5)'~^ + i n—m 



i=2 



(4.39) 



3(5)"-'"-3 3(5)"-"i-3 



54(72) = 2 4 hJm 



5"-™+3 



54! ml 4 



n[ 

i=2 



X 5 + 3 



Pijn - i) 
aAn — i) J 



n[^ 



5J--^+3 



(4.40) 



Here when n ~ m = 1, the products with lower limit two are defined to be one. 



Proof It is clear that Eqs. fl4.39p and fl4.40p are valid for n = m + 1 since /i4(m + 1) = 
8hl{m)s4{m)[5 + 3P4{m)/a4,{m)] and S4,{m + 1) = 8hl{m)sl{m)[7 + /S^^m) / a^^m)] by Eqs. 
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fl4.29p and (14.301) . respectively. Consider Eq. (14.390 holds for a certain positive integer 
n > m, then 

BJn) 



hi{n + 1) = Sh\{n)si{n) 



5 + 3 



5 + 3 



0^4 (n) 



n—m 

m 

i=l 



5 + 3 



Pijn - i) 
a^i^n — i). 



3{5)'-i+l 



X 



J=2 
n—m 

X JJ [5 + 3 

JT, — m + 1 _ 



(3An — i' 



3(5)"-"-3 . 

2 4 /i4(m 



3(5)" 



"S4(m 



5"~'"+3 



a4(n — z) 



2^ 



5 + 3 



'hi{m) 



3(5)" 



^+1 



S4(m) 



c-n — m + 1 



/34(r^ - j) 
ai{n-j) 

n—m 



+3 



i=2 



/34(n - i) 
q;4(?2 — i) 



ai{n — 1) 



5 + 3 



ni-+ 



i=2 



/^4(r^ - j) 
04(72 - j) 



5J-1 



7 + 



3(5)' + l 



/^4(ri - 1) 
ai{n — 1) 



3(5)" 



1 3(5)"-""+l+l gn-m + l„ 

'h^^m) 4 S4(m) 4 



n-m+l ^ / ^ .X 



n—m+l 

' n [ 



5 + 3 



1=1 



a4{n + 1 — i) 



i=2 



ai{n + 1 - j) 



(4.41) 



such that Eq. (14.391) is proved by induction. Eq. (I4.40p can be established by the same 
procedure. □ 

From above lemmas, we have the following bounds for the entropy. 

Lemma IV. 9 The entropy for the number of dimer coverings on SGi{n) is hounded: 
leAmf 



640(5)™[1 - 



f 3 In /i4(m) + In S4(m) + 6 ln2 £4(771) 
< ^SGi — I 



16 



10(5)' 



40(5)' 



■}<0, (4.42) 



where m is a positive integer and e^iji) is defined as 1 — I3i{n) / ai{n) . 



Proof We substitute N{G) = s^in) in Eq. ([2+]) for SG4 so that 

In S4(n) 



S. 



SG4 



lim 



oo5(5" + l)/2 



(4.43) 
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By Lemma flV.St we have 

, , , 3(5)"-™-3, ^ , , 5"-™ + 3, , , 3(5)' 
ms4[n) = lnn.4[m) -\ [ns4[m)-\ 



In 2 + AJn, m) 



(4.44) 



where 

A4(n, m) 



2^ zL ^ 



5 + 3 



i=2 



P^jn - i) 
0:4(72 — i) 



5.-1 + 3 



7 + 



/^4(^ - j) 
04(72 -j) 



(4.45) 



which is bounded as follows. 

By Lemma [IV. 7t we know 64(77.) decreases monotonically to zero for positive integer n. 
€4(1) = 1 - /?4(l)/a4(l) = 4/49. It is easy to find, by Eqs. K36\f and fOTj) . that 

2eAn) 12 



64(11 + 1) 



L8-3e4(n)J 

Therefore, for any integer n > m with m fixed, we have 

6a (771) '^^ 

g,(n + m) = (l + o(n)) . 



Now A4(r;,,m) in Eq. 04.451) can be rewritten 

3[(5)-i - 1] 



A4(r7, m) = 



i=2 



, / 364(^ — 7), 

ln8 + ln(l i^- ^) 

o 



n— m 



E 



3[(5)*-i - 1] r, o 3e4(n-z) 



1=2 

71— m 



ln8 - 



1 + 



3^n,ie4(^ - j) 



(5)^-^ + 3 f 64(n - j) r C/4(n - j) i ^ 

^ 4 1 8 L 8 JJ 

i=i 



where ^n,j,^n,j € (0, 1). It is easy to see that 



A4(n, m) — 



(3 In 2)5" 



5"-'"e4(m),^ , 
-(1 + o{n)) < 777^(1 + o{n)) . 



16 



By Eqs. f OSj) and fOTj) . we have 



^ , , (31n2)5"-"^,^ 

A4(n, m) - ^ { (1 + o{n)) 



(4.46) 



(4.47) 



(4.48) 



(4.49) 



Dimer coverings on the Sierpinski gasket 26 



>_5»-."(i + „(„))(i + ?!lM)"f;'(!lM)^ 

j=0 

e4(m)5"-™(l + o{n)) (l + 
~ 16 — e4(m) 

The proof is completed by taking the infinite n limit in Eq. (14.431) . □ 

The difference between the upper and lower bounds for Ssg4 quickly converges to zero 
as m increases, and we have the following proposition. 

Proposition IV. 2 The entropy for the number of dimer coverings on the four- dimensional 
Sierpinski gasket SG^ln) in the large n limit is Ssd = 0.56337479920.... 

The numerical value of Ssd can be calculated with more than a hundred significant 
figures accurate when m in Eq. fl4.42p is equal to six. It is too lengthy to be included here 
and is available from the authors on request. 



(4.50) 



C. SG^in) 

For the five-dimensional Sierpinski gasket SG^^n), we use the following definitions. 

Definition IV. 3 Consider the five- dimensional Sierpinski gasket SG^ln) at stage n. (i) 
Define f^{n) as the number of dimer coverings such that the six outmost vertices are vacant, 
(a) Define g^{n) as the number of dimer coverings such that one certain outmost vertex is 
occupied by a dimer and the other five outmost vertices are vacant. (Hi) Define hr,[n) as the 
number of dimer coverings such that two certain outmost vertices are occupied by dimers 
and the other four outmost vertices are vacant, (iv) Define r^{n) as the number of dimer 
coverings such that three certain outmost vertices are occupied by dimers and the other three 
outmost vertices are vacant, (v) Define s^{n) as the number of dimer coverings such that 
two certain outmost vertices are vacant and the other four outmost vertices are occupied by 
dimers. (vi) Define t^{n) as the number of dimer coverings such that one certain outmost 
vertex is vacant and the other five outmost vertices are occupied by dimers. (vii) Define 
u^{n) as the number of dimer coverings such that all six outmost vertices are occupied by 
dimers. 
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The quantities f^^n), g^iji), h^{n), r^{n), s^{n), t^{n) and u^{n) are illustrated in Fig. 
[20l where only the outmost vertices are shown. The initial values are /5(0) = 1, (75(0) = 0, 
h^{Q) = 1, r5(0) = 0, 55(0) = 3, ^5(0) = 0, ^5(0) = 15. For the five-dimensional Sierpinski 
gasket SG^ln), the number of vertices is equal to six for n = and odd for all positive 
integer n by Eq. (12.31) . Therefore, h{n), h^{n), s^{n), u^in) are always zero for positive 
integer n. There are (^) = 6 equivalent g^{n) and t^{n), and (3) = 20 equivalent r^{n). 




f^{n) g5{n) h^{n) r^{n) 




S5{n) t^{n) U5{n) 



FIG. 20: Illustration for the dimer coverings f^^n), g^in), h^{n), r5(n), 35(72), t^{n), u^{n). Only the six 
outmost vertices are shown explicitly, where each open circle is vacant and each solid circle is occupied by 
a dimer. 

We write a computer program to obtain following recursion relations. 

Lemma IV. 10 For any positive integer n, 

g^in + 1) = 40gl{n)rr,{n)tl{n) + 5Q0gl{n)rl{n)t5{n) + 424g5{n)rl{n) , (4.51) 

r^in + 1) = 4gl{n)tl{n) + 252gl{n)rl{n)tl{n) + 636g5{n)rl{n)t5{n) + 132rl{n) , (4.52) 
t^in + 1) = A0gl{n)r5{n)tl{n) + 560g5{n)rl{n)tl{n) + 424rl{n)t5{n) , (4.53) 
and for n = 0, 

g,{l) = 280h{0)hl{0)sl{0) + 40f,{0)hl{0)s,{0)u,{0) + 680/1^(0)5^(0) + 2Ahl{0)u,{0) 

= 15840 , (4.54) 

r,{l) = 72f,{0)hl{0)sl{0)u,{0) + 120f,{0)h{0)sl{0) + 712hl{0)sl{0) + 120ht{0)s,{0)u,{0) 
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= 44064 , (4.55) 

t,{l) = 40/5(0)/i5(0)s^(0)M5(0) + 280hl{0)sli0)u5i0) + 24/5(0)5^(0) + 680/i^(0)s^(0) 

= 114912 . (4.56) 

The values of g4{n), r4(n), ^4(72) for small positive n can be evaluated recursively by Eqs. 
f l4.5ip - p:.53p . but they grow exponentially, and do not have simple integer factorizations. 
To estimate the value of entropy for SG^, we define the ratio 

a,{n) = ^ , (4.57) 

and its limit 

as = lim a^i^ri) . (4.58) 

n— >oo 

Lemma IV. 11 Sequence {a5(n)}^j^ decreases monotonically. The limit is equal to 
7399/55. 

Proof From Eqs. (I4.5ip and (14.530 . the ratio t^{n) / g^{n) is invariant. Defined the ratio 
as c, then 

Eqs. (I4.5ip and (14.520 can be modified to be 

g^{n + 1) = %gl{n)r^{n)P^{n) , (4.60) 
r^{n + l)=Agl{n)Q^{n) , (4.61) 

where 

p^[n) = 5c^ + 70cal{n) + 53a^{n) , (4.62) 
Q^(n) = c^ + 63c^al{n) + 159ca^(n) + 33a^(n) . (4.63) 
It is clear that 0:5(72) is bounded below by ^/c because 

rl{n + 1) - g^in + l^in + 1) = rl{n + 1) - cgl{n + 1) 
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= 16[rl{n) - cgl{n)][c'gl{n) + 28c'gl{n)rl{n) + 1542c^gt{n)rl{n) 

+U36cgl{n)rl{n) + 1089r^(n)] > (4.64) 

by induction. It follows that 0:5(77.) decreases for positive n because 

rrojn) _ r^jn + 1) 
g^in) g5{n + l) 

4 



95{n+ 1) 



19rl{n)[rl{n) - cgl{n)] + 53rl{n)[rl{n) - c'gl{n)] + rl{n) - c'gl{n) \ > 

(4.65) 



which implies that the limit 05 exists. From Eqs. (14.601) and (14.611) . we have 



r5(n + 1) + 63c^al{n) + 159caj{n) + 33a 



n] 



(4.66) 



g^{n + l) 2a5(n)[5c2 + 70cai(n) + 53a^(n)] 

By taking the large n limit in Eq. (I4.66P and the requirement that 05 must be real and 
positive, 05 is solved to be ^/c . □ 

The general expressions for gr,{n) and r^{n) in terms of quantities at stage m < n can be 
written as follows. 

Lemma IV. 12 For a non-negative integer m and any positive integer n > m, 

8(6)""^^— 7 — ( — gn — m + 1 j^(^_-^'^n — m gri — m _(^_-i-^'ri — m 

g^{n) = 2 13 g^[rn) 7 r^lm) r 

X H P,{n - z)"-^ n ^^(^ - j)"^^^ , (4.67) 

1=1 j=2 



r5(n) = 2 7 ^5(m) f r5(m) ^ 

n—m . n—m . , 

P,{n - i)^^ \{Q,{n- j) -7 . (4.68) 

i=2 j=l 

Here when n ~ m = 1, the products with lower limit two are defined to be one. 
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Proof It is clear that Eqs. (14.671) and fl4.68p are valid for n = m + 1 since g^im + 1) = 
8(y'|(m)r5(m)P5(m) and r^^m + 1) = 4:g^{m)Q5{m) by Eqs. (14.601) and (14.611) . respectively. 
Consider Eq. (I4.67p holds for a certain positive integer n > m, then 

g^in + l) = 8gl{n)r5{n)P5{n) 

40(6)"-"''-35-5(-l)''-'" ^ 

= 8x2 14 gbim) 7 r^[m)' 

n—m . n—m 

n . 5(6)'-5(-l)' yj 5(6H-l+5(-ip 

i=l j=2 

4(6)"~"'-7+3(-l)"~'" 6"-'"+l-6(-l)"-"' 6"-'"+6(-l)"-"' 

n—m . . n—m 

n6'+6(-l)' -I— r 6J~l-6(-l)^ 
P,{n - t)^^ llQ.in- j) f^P.in) 

i=2 j=l 

= 2 14 gb{f^) ^ ^ 

n—m ^^j^ ^ n—m 

X J] P5(ri - i)'-^^^P,{n - lfP,{n) J] Qsln - j)^^^Q5(^ - 1) 

i=2 j=2 

g^gjn — m-|-l — m + 1 gn — m+2 _j_^_]^~)TT' — "TI + I gn — m + 1 — + 1 

= 2 14 g5{fn) ^ fbi^n) ^ 

n— m+l . n— m+1 . , 

-,— r 6'-(-i)» -I— r 6^-i+(-i)-' 

X n ^5(^+1-^)^^ n Q5(^+i-j) — ^ — , (4.69) 

i=l j=2 

such that Eq. ( 14.67P is proved by induction. Eq. ( 14.68^ can be established by the same 
procedure. □ 

From above lemmas, we have the following bounds for the entropy. 



Lemma IV. 13 The entropy for the number of dimer coverings on SG^ln) is bounded: 

f2lng5{m) lnr5(m) 141n2 Inc 9e5(m) -1 



< 



7(6)"^ 21(6)"^ 21(6)"^ 7(6)™ 56v^(6)' 
279e5(m)2 



448c(6)^ 



1 _ <^5(m) 



(4.70) 



where m is a positive integer and e^{n) is defined as a^{n) — ^/c. 
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A5(n, m) = J2 y In P,{n-z) + J2 / ' In Q,{n - j) , (4.73) 



Proof As 'U5(n) is exactly zero for all positive n, we will substitute N{G) = (75 (n) in 
Eq. for 5G5 so that 

ln^5(n) 
^^^^ = i^3(6^- 

By Lemma [IV. 121 we have 

Qn—m+l _j_ / -^\n—m Qn—m / -^\n—m 

In 5(5 (n) = In 5(5 (m)H lnr5(m) 

+^ In2 + A5(n,m) , (4.72) 

where 

" 6- - (-1)- , "^6^-i + (-l)^- 
In P5(n - z) + 2^ 

i=l j=2 

which is bounded as follows. 

By Lemma [IV.IH we know 65(72) decreases monotonically to zero for positive integer n. 
£5(1) = ^5(1) -^/c= 153/55 - v/399/55. It is easy to find, by Eq. flOTD . that 

e5(n+l) = ^(l + o(n)). (4.74) 

Therefore, for any integer n > m with m fixed, we have 

e,{m + n) = -^^^^{1 + o{n)) . (4.75) 

From Eqs. (14.621) and (14.631) . we have 

P^in + m) = 128c2 + 352cy/ce5{n + m)(l + o{n)) , (4.76) 

and 

Q^{n + m) = 256c^ + dGOc^y/Ze^in + m)(l + o{n)) , (4.77) 

such that 

lnP5(n + m) = 71n2 + 21nc+(l + o(n))ln(^l + i^^^|^^^) , (4.78) 



4v/H 
{n + 



InQ^in + m) = 8 ln2 + 3 Inc + (1 + o(n)) In ( 1 + . (4.79) 

V 4a/c / 

Now A5(n,m) in Eq. (14.731) can be rewritten, 

A5(n,m) = ^Ml=izll![7in2 + 21nc] 

1=1 
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i=l 

n—m 



J=2 



(4.80) 



It is easy to see that 



A.(«.m)- [H|i^ + i^]6"-»(l +o(„))>HI|M6..-".(l -<,(„)) , (4,81) 

By Eqs. f lOHD and fHTHl) . we have 

A / N rl01n2 31nc] „ ^, , 
A5(n, m) - + 6"-™(l + o(n)) 

i=o 

27e5(m)6"-"^(l + o(n)) (^1 + i^gg^^ 
- 7[8v/5-e5(m)] ^ " ^^"^^^ 

The proof is completed by taking the infinite n hmit in Eq. (14.71 1) . □ 

The difference between the upper and lower bounds for Ssgs quickly converges to zero 
as m increases, and we have the following proposition. 



Proposition IV. 3 The entropy for the number of dimer coverings on the five- dimensional 
Sierpinski gasket SG^i^n) in the large n limit is Ssgs = 0.67042810305.... 

The numerical value of Ssgs can be calculated with more than a hundred significant 
figures accurate when m in Eq. fl4.70p is equal to six. It is too lengthy to be included here 
and is available from the authors on request. 

We notice that the convergence of the upper and lower bounds of the entropy for dimer 



coverings on SGd{n) is about the same for = 3, 4, 5, similar to the results observed in [27 1 
for the dimer-monomer model on SGd{n). 
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V. SUMMARY 



Compare the present results with those in Ref. [27|], it is clear that the number of 
dimer coverings on the Sierpinski gasket SGd{n) is less than that of dimer- monomers. The 
asymptotic growth constant zscat ^'^^ dimer-monomer model defined as Eq. (2.1) of 



27l | corresponds to the entropy Ssc^t ^'^^ dimer coverings defined in Eq. (12. ip . We 
summarize the values of SsCdb ratio SsGab/^SGat Table [H The value of Ssg^ 

increases as dimension d increases. Similarly for the generalized two-dimensional Sierpinski 
gasket, the exact value of 5*5^2 increases slightly as b increases. For the cases studied, the 
ratio SsGdl^SGd also increases as dimension d increases, and SsG2b/^SG2b increases slightly 
as b increases. 

It is interesting to compare entropy of dimer coverings on the Sierpinski gasket SG^ with 
that on the (i-dimensional hypercubic lattice which is also 2(i-regular. The entropy of 
the square lattice was known to be G/tt ^, where G is the Catalan number, for decades, 
while the entropy of the simple cubic lattice was estimated to be 0.44647 [28]. They are 
relatively larger than the entropies on SGd with d = 2,3 presented here. The values of Sc^^ 
and the ratio SsGa/^Cd d = 2,3 are given in Table [H It appears that as the d increases, 
the value Ssg^ approaches to the value Sc^ from below. As we have obtained the highly 
accurate value for the entropy on SG^ with d = 4, 5, there is no numerical estimation for 
the entropy on C^. with d > 4, to the best of our knowledge. 
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